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The  purpose  of  the  current  paper  is  to  investigate  the  Cynara  cardunculus  (cardoon)/lignite  co-firing 
concept  at  a  10%  biomass  thermal  share  for  a  330  MWe  pulverized  fuel  plant  in  Northern  Greece. 
Two  main  cases  are  considered:  (A)  co-milling  of  pelletized  biomass  with  lignite  and  (B)  milling  of 
balled  and/or  shredded  material  in  dedicated  biomass  mills.  A  CFD  model  for  plant  simulation  is 
developed  taking  into  account  the  particularities  of  co-firing  operation,  in  particular  the  non-spherical 
form  of  the  biomass  particle  and  its  effect  on  the  drag  coefficient,  devolatilization  and  combustion 
rates.  Simulations  are  performed  for  the  reference  case  of  lignite-firing,  which  show  good  agreement 
with  plant  operational  data,  and  for  the  two  co-firing  cases.  Further  subdivision  of  the  co-firing  cases 
is  based  on  the  biomass  particle  size  and  injection  point  in  the  furnace.  In  most  co-firing  cases,  CFD 
results  indicate  that  the  substitution  of  lignite  with  biomass  has  minimal  impact  on  the  plant  oper¬ 
ational  parameters  as  well  as  the  potential  for  NOx  reductions.  Increased  unburnt  losses  in  the  fly  and 
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Nomenclature 

A 

pre-exponential  factor  in  Arrhenius  expression  for  char 

n 

refractive  index,  dimensionless 

combustion,  kg  m”2  s”1  Pa-1 

klreac 

reaction  order,  dimensionless 

^bio 

pre-exponential  factor  for  devolatilization  mechanism 

Nu 

Nusselt  number,  Nu  =  h  dp\k00 

of  biomass,  s”1 

NPV 

Net  Present  Value 

Alig 

pre-exponential  factor  for  devolatilization  mechanism 

Op 

heat  exchanged  between  the  particle  and  the  gas  phase, 

of  lignite,  s”1 

j 

Ap 

heat  transfer  surface  area  of  the  particle,  m2 

r 

discount  rate,  % 

^sph 

surface  area  of  the  droplet,  m”2 

R 

kinetic  rate,  kg  m"2  s”1  Pa”1 

Bi 

Biot  number,  dimensionless 

Rg 

ideal  gas  constant,  J  kmol”1  K”1 

C 

linear-anisotropic  phase  function  coefficient,  dimen- 

Re 

relative  Reynolds  number,  dimensionless 

sionless 

t 

time,  s  or  year 

CD 

drag  coefficient,  dimensionless 

T 

gas  temperature,  K 

CFt 

net  cash  flow  of  the  investment  within  a  period  t 

T 

lp 

particle  temperature,  I< 

9 

specific  heat  per  unit  mass,  kj  kg-1  K”1 

TPC 

Total  Plant  Cost 

D0 

mass  diffusion  limited  rate  on  char  burnout,  kgm”2 

T 

1  00 

temperature  of  the  continuous  phase,  I< 

s”1  Pa”1 

Ui 

fluid  phase  velocity  for  the  i  direction  in  Cartesian  coor- 

dp 

particle  diameter,  m 

dinates,  m  s”1 

E 

activation  energy  for  char  combustion,  J  kmol”1 

Up 

particle  phase  velocity,  m  s”1 

Ebio 

activation  energy  for  devolatilization  mechanism  of  bio¬ 
mass,  J  kmol”1 

v 

1  ox 

local  mass  fraction  of  oxidant  in  the  gas 

Eiig 

activation  energy  for  devolatilization  mechanism  of  lig- 

Greek  letters 

nite,  J  kmol”1 

a 

absorption  coefficient,  m”1 

fh 

fraction  of  heat  from  char  combustion  absorbed  by  the 

8p 

particle  emissivity,  dimensionless 

particle,  dimensionless 

0 

shape  factor,  dimensionless 

G 

incident  radiation,  W  m”2  sr”1 

radiation  temperature,  K 

gi 

gravitational  acceleration  for  the  i  direction  in  Cartesian 

P 

molecular  viscosity  of  the  fluid,  kg  m”1  s”1 

coordinates,  m  s”2 

convective  heat  transfer  coefficient,  W  m”2  K”1 

P 

fluid  phase  density,  kg  m”3 

h 

Pp 

density  of  particle,  kg  m”3 

hfg 

latent  heat  of  vapourization,  J/kg 

G 

Stefan  Boltzmann  constant,  5.67  x  10”8  W  m”2  K”4 

Hreac 

heat  released  by  the  char  combustion  reaction,  J/kg 

scattering  coefficient,  m”1 

IRR 

Internal  Rate  of  Return,  % 

1 fieri 

enhancement  factor,  dimensionless 

mp 

mass  of  the  particle,  kg 

Mw,ox 

molecular  weight  of  oxides,  kg  kmol”1 

N 

duration  of  investment,  year 

bottom  ash  can  be  expected  when  large  biomass  particles  are  combusted;  the  best  scenario  appears 
to  be  the  separate  milling  of  biomass  and  its  injection  in  the  vapour  burners.  Economic  evaluation  of 
the  two  cases  also  favours  the  dedicated  milling  system,  despite  its  higher  investment  cost,  due  to  the 
lower  fuel  cost  of  non-pelletized  biomass.  The  co-firing  concept  appears  to  be  highly  profitable  for  the 
plant  operator  and  further  optimisation  of  the  combustion  process  can  be  suggested  by  the  installa¬ 
tion  of  dedicated  biomass  burners. 

©  2013  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

The  improvement  of  the  efficiency  and  the  environmental  per¬ 
formance,  including  reduction  of  greenhouse  gas  emissions,  of  fos¬ 
sil  fuel  utilisation  for  energy  production  has  led  to  significant  R&D 
efforts  in  the  field  of  novel  combustion  technologies,  such  as  oxy- 
fuel  combustion  [1]  and  flameless  oxidation  furnaces  [2].  However, 
most  of  these  technologies  are  in  early  stages  of  development  and 
their  implementation  in  existing  fossil  fuel  power  plants  faces  sev¬ 
eral  limitations.  Part  of  the  solution  in  the  quest  for  reducing  GHG 
emissions  from  existing  coal-fired  power  plants  can  be  offered  by 
biomass  co-firing,  a  low  cost  technology  option  [3],  the  importance 
of  which  in  reaching  the  EU  targets  for  the  increase  of  the  share  of 
renewable  energy  sources  in  the  energy  sector  is  well  recognised 
by  EU  authorities  [4]  and  research  academy  [5].  The  most  common 
technological  option  is  direct  combustion  of  the  solid  biofuels, 
although  examples  of  syngas  combustion  after  biomass  gasifica¬ 
tion  are  also  presented  [6]. 

Co-firing  is  a  particularly  promising  solution  for  GHG  emissions 
decrease  in  Greece,  due  to  the  high  installed  capacity  of  lignite-fired 


power  plants  and  the  large  contribution  of  indigenous  lignite  in  the 
gross  energy  production  [7].  However,  instead  of  the  usual,  typi¬ 
cally  high  quality  biogenic  feed  stocks  used  in  most  of  the  co-firing 
power  plants  across  Europe  (saw  dust  or  wood  pellets),  Greece  will 
have  to  rely  mostly  on  the  locally  produced  agricultural  biomass  for 
reasons  of  decreasing  the  cost  of  the  supply  chain  [8].  Despite  the 
high  availability  of  agricultural  residues  in  Greece  [9],  the  lack  of 
mechanization  for  residue  harvesting,  the  associated  high  operating 
costs  and  other  competing  uses  have  shifted  attention  towards  the 
exploitation  of  locally  grown  energy  crops  for  power  generation. 
One  of  these  crops  to  which  Greek  stakeholders  have  paid  attention 
during  the  last  years  is  Cynara  caxdunculus  L.  (cardoon),  a  perennial, 
herbaceous  crop  of  Mediterranean  origin  well  adapted  to  the 
xerothermic  conditions  of  Southern  Europe  [10].  This  is  owed  to 
its  potential  for  high  yields  [11].  Cardoon  is  currently  cultivated 
in  the  Region  of  Western  Macedonia  in  the  frame  of  a  regional 
development  program,  aiming  to  produce  biomass  for  combustion 
reasons  for  the  nearby  Kardia  power  plant  [12]. 

An  overview  of  co-firing  experiences  have  been  collected  in  re¬ 
view  papers  by  [13,14].  Several  studies  concerning  the  modelling 
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of  co-firing  solid  biomass  along  with  coal  in  pilot  reactors  or  large- 
scale  pulverized  boilers  have  been  published  [15-21].  A  key  factor 
outlined  in  most  of  the  simulations  is  the  large  length  to  diameter 
ratio  of  biomass  particles,  which  influences  their  trajectories  inside 
the  boiler  and  consequently  their  residence  time  and  final  burnout 
[16].  The  co-firing  of  cardoon  with  a  hard  coal  blend  was  investi¬ 
gated  by  [19]  for  a  front-fired  industrial  furnace  with  an  output 
of  350  MWe,  showing  that  is  a  promising  fuel. 

The  purpose  of  the  present  study  is  to  evaluate  the  effect  of  co¬ 
firing  cardoon  with  Greek  lignite  in  a  330  MWe  super-critical,  Ben¬ 
son  type,  tangentially  fired  boiler,  taking  into  account  the  major 
physical  and  chemical  differences  in  the  combustion  process  for 
cardoon  and  lignite.  The  optimal  firing  scheme  is  the  main  topic 
investigated  in  this  paper,  intending  to  add  value  to  the  investiga¬ 
tions  related  to  cardoon  specifically  or  to  biomass  as  an  alternative 
fuel  in  general.  The  main  parameters  to  be  investigated  are  related 
to  fuel  burnout,  temperature,  gas  species  and  NO*  spatial  distribu¬ 
tions  for  different  firing  concepts.  The  numerical  results  are  evalu¬ 
ated  on  the  basis  of  the  aforementioned  values  at  the  furnace  exit, 
within  the  furnace  and  at  the  hopper  outlet.  Numerical  modelling 
of  tangentially  fired  boilers  have  been  performed  by  several 
authors  [22-27],  including  studies  on  furnaces  operating  with 
low  quality  coals  [28-30].  However,  co-firing  studies  with  lignite 
on  tangentially-fired  boilers  have  been  limited  to  solid-recovered 
fuels  [31]. 

In  addition,  an  economic  evaluation  of  the  investigated  co-firing 
is  performed  by  taking  into  account  all  relevant  project  costs  (Cap¬ 
ital  Expenditure  -  CAPEX,  Operational  Expenditure  -  OPEX)  and 
calculating  the  main  economic  indices  (Net  Present  Value  -  NPV, 
Internal  Rate  of  Return  -  IRR,  payback)  for  different  cost  cases  of 
biomass  and  Green  House  Gas  emission  allowances.  Combined, 
the  two  approaches  offer  an  overall  techno-economic  optimisation 
of  the  biomass  feeding  and  firing  system  for  the  investigated  power 
plant. 

2.  Methodology 

Following  the  approach  presented  in  the  Introduction,  the 
methodology  section  is  divided  into  two  sections.  In  the  first  part, 
the  methodology  for  the  technical  evaluation,  i.e.  the  numerical 
simulations  of  the  industrial  boiler  (CFD),  is  presented.  More  spe¬ 
cifically,  the  basic  numerical  models  applied  in  the  CFD  code  are 
presented  and  the  specific  models  developed  for  modelling  the 
combustion  behaviour  of  biomass  particles  are  thoroughly  de¬ 
scribed.  In  the  second  section,  the  methodology  for  the  economic 
evaluation  is  presented. 

2.2.  CFD  methodology 

2.2 .2.  Basic  numerical  models 

The  computational  domain  includes  the  main  combustion  zone 
of  the  furnace;  it  excludes  the  convective  part  of  the  boiler.  The 
implemented  numerical  mesh  is  built  using  the  GAMBIT  software 
and  consists  of  940.315  unstructured  tetrahedral  cells.  The  grid  is 
locally  refined  in  the  near-burner  area,  the  hopper  region  and  the 
flue  gas  suction  ducts  in  order  to  increase  the  accuracy  of  the 
numerical  results  near  the  boundary  walls,  since  the  heat  flux 
transferred  to  the  walls  is  a  crucial  parameter  for  the  evaluation 
of  the  boiler  operation.  The  grid  density  is  comparable  with  the 
meshes  employed  in  similar  simulations  of  full-scale  boilers 
[1,19,24,32,33].  This  grid  density  has  as  also  been  validated  as 
dense  enough  by  performing  simulations  with  a  coarser  numerical 
grid  of  around  700.000  cells;  the  deviations  in  the  calculated  flue 
gas  characteristics  at  the  furnace  outlet  are  around  0.35%  com¬ 
pared  to  the  ones  presented  in  this  work.  The  smallest  cell  volume 
is  equal  to  4.42E-05  m3  and  the  largest  one  equal  to  0.268  m3. 


The  numerical  simulations  are  performed  using  the  commercial 
CFD  code  Fluent  v.12.1,  [34]  in  a  Pentium  4  2.4  GHz  personal  com¬ 
puter.  Each  run  requires  about  3  days  of  real  computing  time  with 
parallel  processing  (4-CPU  cores),  including  the  NOx  simulations. 
The  standard  time  averaged  mathematical  expressions  of  equa¬ 
tions  for  mass,  momentum,  enthalpy  and  species  are  solved  using 
the  SIMPLE  algorithm  for  predicting  the  flow,  temperature  and 
concentration  of  gas  species  within  the  boiler.  All  variables  (veloc¬ 
ity  components,  temperature,  and  gas  species  including  NOx)  are 
solved  using  a  second  order  Upwind  scheme.  The  standard  k-e 
model  is  used  along  with  standard  wall  functions  for  the  flow  char¬ 
acteristics  near  the  boundary  conditions  of  boiler’s  wall. 

Normally,  the  values  of  wall  emissivity  and  wall  temperature 
exhibit  a  spatial  as  well  as  a  temporal  distribution  because  of  the 
different  deposition  rates  of  ash  particles  onto  the  heat  exchangers. 
For  a  more  valid  simulation  of  the  phenomenon,  one  should  mea¬ 
sure  these  rates  either  by  pyrometers  or  thermal  imaging  cameras 
and  use  the  extrapolated  data  as  an  input  in  the  numerical  model. 
In  the  present  study,  due  to  the  lack  of  such  data  for  this  large  scale 
industrial  boiler,  standard  values  are  adopted.  The  wall  tempera¬ 
ture  is  set  to  a  temperature  equal  to  365  °C,  i.e.  50°  higher  than 
the  average  steam  temperature  in  the  steam  tubes  [19],  whilst  a 
typical  value  of  0.7  for  wall  emissivity  is  used  [1,25,30,35-38]. 
The  wall  emissivity  value  is  found  to  produce  accurate  results  of 
the  wall  heat  flux  and  furnace  exit  temperature  for  the  reference 
case  of  lignite  combustion. 

The  PI  radiation  model  for  the  simulation  of  the  radiation  heat 
transfer  is  used  and  the  absorption  coefficients  of  the  gas  phase  are 
calculated  using  the  domain  based  weighted-sum-of-grey-gases 
model  (WSGGM),  [34].  The  transport  equation  for  the  incident 
radiation  is  given  below: 

V  •  (rVG)  -  a  •  G  +  4  •  a  •  n2  •  cr  •  T4  =  0,  (1) 

where  the  parameter  r  is  given  by  the  equation: 


(3  ■  (a  +  os)  -  C  ■  gs) 

The  impact  of  the  discrete  phase  particles  with  the  radiation 
model  has  not  been  considered  in  Eq.  (1)  in  this  work.  The  effect 
of  the  lignite  particles  on  the  radiation  model  is  expected  to  be 
negligible,  since  both  the  particle  emissivity  and  the  convective 
heat  transfer  coefficient  experienced  by  the  particles  are  quite 
high:  0.9  and  around  4500  W/m2  K  respectively.  As  a  result,  the 
temperature  of  the  lignite  particles  is  expected  to  be  close  to  that 
of  the  surrounding  gas.  The  larger  biomass  particles  have  a  poten¬ 
tially  larger  impact  on  the  incident  radiation  transport  equation; 
therefore,  these  terms  should  be  considered  in  future  works,  espe¬ 
cially  for  higher  biomass  thermal  shares  or  dedicated  biomass 
combustion.  In  the  present  study  however,  the  deviation  is  mini¬ 
mal  due  to  the  low  share  of  biomass  particles  in  the  overall  particle 
load  of  the  boiler. 

The  furnace  exit  and  the  flue  gas  recirculation  ducts  are  mod¬ 
elled  using  outflow  boundary  conditions.  The  relative  distribution 
of  flue  gases  between  the  recirculation  ducts  and  the  furnace  exit 
is  based  on  the  heat  and  mass  balance  calculations  for  the  milling 
system.  Typically,  16.6-17.6%  of  the  flue  gas  leaves  the  furnace 
through  the  flue  gas  recirculation  ducts. 

An  Eulerian-Lagrangian  description  is  adopted  to  resolve  the 
coupling  of  gas  with  the  particulate  phase.  The  equations  for  the 
description  of  the  two  discrete  phases  (lignite  and  cardoon  parti¬ 
cles)  are  described  in  detail  in  the  following  paragraphs.  Homoge¬ 
neous  gas  phase  reactions  between  species  released  during  fuel 
devolatilization  are  modelled  using  the  Finite  Rate/Eddy  Dissipa¬ 
tion  concept.  Radiation  and  chemistry  are  updated  every  one  (1) 
gas  flow  iteration.  In  the  frame  of  this  work,  the  reactions  and 
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the  corresponding  Arrhenius  rate  coefficients  used,  are  based  on  a 
Westbrook  and  Dryer  2-step  global  combustion  mechanism  for 
hydrocarbon  fuels  [16,20,39].  The  mechanism  is  applied  to  two 
“artificial”  species,  LV  and  BV,  representing  the  lignite  and  biomass 
volatiles  respectively.  Table  1  summarises  the  relevant  reactions 
and  associated  rates. 


2.2.2.  Discrete  phase  modelling 

Fuel  particles  are  simulated  in  a  Lagrangian  fashion  and  particle 
trajectories  are  calculated  throughout  the  whole  computational 
domain.  Particle  dispersion  by  the  induced  turbulence  of  the  gas 
phase  is  simulated  using  the  stochastic  tracking  model,  which  in¬ 
cludes  the  effect  of  instantaneous  turbulent  velocity  fluctuations 
of  the  gas  on  the  particle  trajectories.  The  flow/particle  interactions 
are  represented  by  appropriate  sink/source  terms  both  in  the  gas/ 
particulate  equations,  while  these  terms  are  updated  once  every  25 
iterations  of  the  gas  fluid  flow. 

In  co-firing  simulations,  the  different  behaviour  of  lignite  and 
biomass  particles  deserves  special  attention.  In  the  CFD  simula¬ 
tions  coal  and  lignite  particles  are  considered  as  spherical  and 
the  equations  of  motion  and  combustion  modes  they  obey  can  be 
handled  by  the  pre-built  models  of  most  of  the  commercial  CFD 
packages.  On  the  other  hand,  biomass  particles  are  highly  irregular, 
and  their  approximation  as  spherical  increases  the  deviation  from 
real  conditions.  To  reduce  this  deviation,  in  most  cases  they  can  be 
regarded  as  having  a  cylindrical  shape.  Thus,  this  affects  both  the 
induced  drag  force  by  the  gas  phase,  while  also  increasing  the  heat¬ 
ing,  devolatilization  and  combustion  rates.  A  parameter  depicting 
the  irregularity  of  solid  particles  is  the  shape  factor  <9,  which  is 
the  particle  surface  area  of  a  sphere  having  the  same  volume  as 
the  particle  divided  by  the  actual  surface  area  of  the  particle.  For 
spherical  lignite  particles,  it  is  obvious  that  0  is  equal  to  1,  while 
for  biomass  particles  a  constant  value  of  ©  equal  to  0.78  is  used 
in  the  present  study.  This  value  is  derived  for  a  cylindrical  shaped 
biomass  particle,  which  length  to  diameter  ratio  is  equal  to  4.  The 
value  of  4  is  a  quite  often  value  found  in  the  literature  and  stands 
also  for  the  investigated  cardoon  particles  [16,19]. 

The  particle  equation  of  motion  is  solved  for  each  trajectory. 
In  order  to  ensure  good  statistics,  a  total  of  102,050  particle- 
trajectories  are  simulated.  The  equation  of  motion  for  each 
particle  is: 


dUp 

dt 
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For  lignite  particles,  the  drag  coefficient  is  given  by  the  standard 
equation  of  [40],  while  for  biomass  particles,  the  effect  of  non¬ 
sphericity  on  the  drag  coefficient  is  modelled  using  the  equations 
of  [41],  which  is  the  default  Fluent  model  for  non-spherical 
particles. 

The  fuel  particles  undergo  five  sequential  heterogeneous  steps: 
inert  heating  till  the  particle  reaches  the  vapourization  tempera¬ 
ture;  drying;  inert  heating  till  the  particle  reaches  the  devolatiliza¬ 
tion  temperature;  devolatilization;  char  burnout  and  a  final  inert 


heat  step  for  the  remaining  ash  particle.  The  calculations  of  the 
devolatilization  and  combustion  rates  of  lignite  and  biomass  parti¬ 
cles  can  be  performed  by  several  sophisticated  models,  such  as  the 
Intrinsic  Model,  which  take  into  account  the  chemical,  molecular 
structure  and  the  intrinsic  surface  of  solid  fuels.  For  this  work, 
the  chemical  rate  is  calculated  based  on  simpler  Arrhenius  models, 
due  to  restrictions  on  the  availability  of  data  for  the  investigated 
fuels. 

For  the  devolatilization  mechanism,  a  single  rate  model  ex¬ 
pressed  in  an  Arrhenius  form  is  adopted.  For  lignite,  the  pre¬ 
exponential  factor  is  equal  to  Aug=  5  x  1013  (s-1)  and  the  activation 
energy  equal  to  Eiig=  7.4  x  107  (J/kmol)  [42].  For  biomass,  the  cor¬ 
responding  rates  are  Abi0  =  9  x  1018(s-1)  and  Ebio  =  (2.39  x  1 0s  J / 
kmol)  [19]. 

Char  combustion  is  modelled  using  a  kinetic/diffusion  limited 
model.  Specifically  for  the  Greek  lignite  particles  the  char  combus¬ 
tion  model  and  the  kinetic  parameters  for  two  independent  paral¬ 
lel  reactions  are  described  by  Vamvuka  et  al.  [43].  This  model  has 
also  been  implemented  in  the  numerical  model  by  using  custom 
built  functions,  so  as  the  results  are  more  representative  for  the 
Greek  lignite  particles.  For  biomass  combustion,  the  kinetics  of 
the  char  combustion  are  modelled  using  a  single  step  reaction 
the  controlling  parameters  being  equal  to  A  =  0.3272  (kg/ 
(m2  s  Pa))  and  E  =  4.5758  x  107  (J/kmol)  [19].  The  overall  biomass 
combustion  rate,  based  on  the  kinetics/diffusion  limited  model  is 
given  by  the  following  modified  equation  [34],  where  the  standard 
expression  is  multiplied  by  the  factor  c/>en : 

dvcip _  .  p  •  Rg  •  T oo  •  Y ox  D0  R. 
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The  burning  enhancement  factor,  </>en,  was  developed  by  Gera 
et  al.  [16]  in  order  to  take  into  account  the  fact  that  there  is  an  in¬ 
crease  in  the  overall  burning  rate  of  a  cylindrical  biomass  particle 
compared  to  a  spherical  one  of  the  same  surface.  The  enhancement 
factor  is  given  by  the  following  simple  algebraic  function: 


0.3  -0  +  0.7 
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The  particle  temperature  Tp  is  calculated  by  considering  the 
heat  transfer  due  to  convection  and  radiation,  as  well  as  the  heat 
exchange  with  the  surrounding  gas  during  each  step  (drying,  dev¬ 
olatilization,  char  combustion)  using  the  following  equation: 


d(Cp  •  Tp) 
dt 
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The  term  Qp  in  Eq.  (4)  refers  to  the  heat  exchanged  between  the  par¬ 
ticle  and  the  gas  phase  through  water  vapourization,  Eq.  (5),  or 
through  the  heterogeneous  char  combustion,  Eq.  (6): 
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Table  1 

Kinetic  rates  for  gas  phase  reactions. 


Reaction 

A 

.  n-1  la 

L  read 

E  (J/kmol) 

Reaction 

order5 

Westbrook  and  Dryer  2-step  global  combustion  mechanism  for  hydrocarbon  fuels 

LV  +  02  -►  CO  +  H20  +  S02  +  N2 

2.119  x  1011 

2.027  x  108 

[LV]02  [02]13 

BV  +  02  CO  +  H20  +  S02  +  N2 

2.119  x  1011 

2.027  x  108 

[BV]°-2-[02]13 

CO  +  0.5  02  ->  CO 

2.239  x  1012 

1.702  x  108 

[CO]  [02]025 
[H20]os 

LV:  lignite  volatiles,  BV:  biomass  volatiles. 
a  nreac:  Total  reaction  order. 
b  Concentration  values  in  (kmol/m3). 


where  hfg  is  the  latent  heat  of  vapourization,  Hreac  is  the  heat  re¬ 
leased  by  char  combustion  and  fh  is  the  heat  fraction  observed  by 
the  particle.  In  this  work,  the  char  burnout  product  is  assumed  to 
be  C02  and  the//,  is  set  to  0.3  [34,44]. 

It  should  be  noted  that  the  effect  of  non-sphericity  for  the  bio¬ 
mass  particles  is  taken  into  account  in  Eq.  (6)  by  using  the  actual 
surface  area  of  the  particle  and  not  that  of  the  spherical  equivalent. 
Moreover,  for  the  cylindrical  particles,  the  heat  transfer  coefficient 
is  evaluated  using  the  well  suited  for  this  case  Churchill  and  Bern¬ 
stein  expression  for  the  Nusselt  number  [45].  Intra-particle  con¬ 
duction  of  heat  and  mass  are  not  considered  in  the  present 
study.  For  the  lignite,  the  calculation  of  the  Biot  number  ( =h  Deq\ 
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k)  yields  values  of  1.06,  1.51  and  1.48  for  particles  diameters  of 
50  pm,  90  pm  and  170  pm  respectively.  Normally,  for  such  values 
conduction  within  the  particle  volume  is  a  limited  factor,  but 
according  to  [46-50]  the  heat  transfer  is  limited  by  the  internal 
conduction  when  Biot  number  is  larger  than  10.  Therefore,  for 
the  lignite  particles  the  neglect  of  heat  conduction  mechanism 
within  the  particle  is  not  a  coarse  approach,  but  for  the  case  of  bio¬ 
mass  particles  this  effect  should  be  considered,  as  in  the  approach 
of  [46],  especially  for  cases  when  the  biomass  thermal  share  in  the 
fuel  mixture  is  high. 

The  char  combustion  model  employed  is  based  on  the  work  of 
Smith  [50],  which  assumes  that  oxygen  fully  penetrates  the  pores 
of  the  particles  so  that  the  particle  diameter  remains  constant  and 
the  density  is  decreasing  during  combustion. 

Due  to  limitations  of  the  standard  sub-models  of  the  ANSYS/ 
Fluent  code,  all  steps  one  biomass  particle  undergoes  during  its 
combustion  are  implemented  in  the  ANSYS  /Fluent  standard  plat¬ 
form  by  the  use  of  custom  built  User  Defined  Functions  (UDFs). 

To  ensure  the  solution  convergence  the  dimensionless  residuals 
for  the  continuity  equation  should  be  less  than  10”4,  while  for  the 
momentum,  energy  and  NOx  calculations  residuals  should  be  less 
than  10-6  after  the  introduction  of  all  submodels.  Moreover,  the 
area  weighted  averaged  values  of  the  main  flue  gas  species  concen¬ 
trations  and  temperature  at  three  representative  horizontal  levels 
of  the  boiler  was  monitored;  if  their  values  remains  constant  with 
less  than  0.3%  variation,  then  the  solution  is  assumed  to  have 
reached  convergence. 

2 A  3.  NOx  formation  modelling 

In  the  present  study,  both  thermal  and  fuel  NO*  mechanisms  are 
considered.  Prompt  NOx  formation  is  not  taken  into  consideration, 
since  it  is  only  relevant  in  fuel-rich  conditions,  uncharacteristic  of 
pulverized  fuel  systems. 

For  thermal  NOx  formation  the  extended  Zeldovich  mechanism 
is  used.  The  mass  transport  equation  for  the  NO,  FICN  and  NH3  spe¬ 
cies  including  their  convection,  diffusion,  production  and  con¬ 
sumption  mechanisms  are  solved.  The  transport  equation  for  NO 
is  solved  by  estimating  the  local  concentrations  of  O  and  OH  radi¬ 
cals  through  the  partial  equilibrium  approach  and  taking  into  ac¬ 
count  the  effect  of  turbulence  interaction  by  the  probability 
density  function  (PDF)  [34]. 

Regarding  the  fuel  NOx  mechanism,  nitrogen  oxides  are  pro¬ 
duced  by  the  oxidation  of  the  nitrogen  bound  in  the  coal,  both  in 
the  volatile  matter  and  in  the  char.  The  HCN  and  NH3  generally  re¬ 
act  with  02  to  form  either  NO  in  fuel-lean  regions  or  with  NO  to 
form  N2  in  fuel-rich  regions  [51]. 

Fuel  bound  N  for  the  two  examined  fuels  (lignite  and  cardoon) 
is  considered  to  be  distributed  between  volatiles  and  char  at  a  ratio 
equal  to  the  respective  ratio  of  volatile  and  char  content  of  the  fuel. 
The  nitrogen  containing  species  released  from  volatiles  are  HCN 
and  NH3,  with  a  ratio  of  9:1  for  lignite  [52]  and  1:9  for  biomass 
[53].  Char  bound  nitrogen  is  considered  to  directly  convert  to 
NO,  mainly  as  a  desorption  product  from  oxidised  char  nitrogen 
atoms  [54].  A  conversion  factor  of  0.7  is  adopted  for  the  char  bound 
nitrogen  [55].  The  NO  formed  by  the  previously  described  mecha¬ 
nisms  can  be  reduced  via  its  surface  reaction  with  char  particles 
[56].  The  detailed  equations  of  the  NOx  formation  and  destruction 
mechanisms  used  in  this  work  can  be  found  in  [1]. 

2.2.  Economic  evaluation  methodology 

For  the  evaluation  and  the  comparison  of  the  two  investigated 
cases  which  represent  different  technical  solutions  for  the  adop¬ 
tion  of  co-firing  at  the  Meliti  PP,  the  Net  Present  Value  of  invest¬ 
ment  (NPV),  the  Internal  Rate  of  Return  (IRR)  and  the  repayment 
period  of  investment  are  used  as  criteria. 


The  Net  Present  Value  (KPA-Net  Present  Value/NPV)  is  the  sum 
of  all  cash  flow  after  taxes  deduced  in  the  reference  year.  It  is  an 
indicative  of  value  of  money  through  time.  Investments  with  posi¬ 
tive  NPV  are  profitable  and  favourable.  The  NPV  is  calculated 
according  to  the  equation: 
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where  NPV  is  the  Net  Present  Value  of  investment,  CFt  is  the  net 
cash  flow  of  the  investment  within  period  t,  r  the  discount  rate 
(%)  and  t  the  time  period,  which  takes  values  from  0  to  N  (years). 
N  is  the  length  of  the  investment. 

The  Internal  Rate  of  Return  (IRR)  is  defined  as  the  money  mar¬ 
ket  rate  at  which  the  NPV  is  zero,  which  means  that  the  present 
value  of  capital  investment  equals  the  present  value  of  revenues 
from  the  operation.  The  IRR  is  calculated  from  the  following 
equation: 


°  =  E 
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The  main  feature  of  IRR  is  that  it  can  be  used  for  investments  of 
different  sizes.  If  the  calculated  IRR  of  a  project  is  higher  than  the 
discount  rate,  then  the  investment  is  considered  as  profitable, 
otherwise  the  investment  is  discarded.  The  higher  the  IRR,  the 
more  profitable  the  project  is. 

The  repayment  period  of  the  investment  is  the  time  required  to 
recover  the  initial  investment  and  can  be  calculated  simply  by 
dividing  the  initial  investment  with  the  net  profit/losses.  In  this 
work,  the  repayment  period  is  calculated  by  taking  into  account 
the  change  in  value  of  money  over  time. 


3.  Case  studies 

Similar  to  the  way  of  approach  in  the  previous  parts,  this 
part  is  divided  in  two  sections.  In  the  first  section  the  operating 
and  boundary  conditions  of  the  CFD  simulations  are  presented 
and  the  different  cases  considered  for  evaluation  are  analysed. 
In  the  second  section  the  detailed  description  of  the  technical 
solution  for  each  considered  case  is  presented,  so  that  an  esti¬ 
mation  of  capital  and  operating  expenditures  for  each  case  is 
possible. 

Generally,  two  main  cases  are  considered  for  both  CFD  model¬ 
ling  and  economic  analysis: 

Case  A  foresees  the  delivery  of  biomass  in  the  lignite  yard  in  the 
form  of  pellets.  The  technical  set-up  includes  facilities  for  inter¬ 
mediate  storage  and  feeding  of  biomass  to  the  main  lignite  con¬ 
veyor  belt.  Following  that,  the  biomass  fuel  is  handled  by  the 
existing  system  for  lignite,  e.g.  it  passes  through  the  same 
crushers  and  mills  and  enters  the  boiler  through  the  existing 
burners  without  any  further  optimisation. 

In  Case  B,  biomass  is  delivered  to  the  power  plant  as  baled 
material  or  in  bulk,  as  shredded  particles.  Such  forms  cannot  be 
handled  by  the  existing  feeding  system  for  lignite  and  thus  a  need 
for  a  dedicated  biomass  handling  line  arises.  The  key  point  for  the 
CFD  simulations  is  that  biomass  is  milled  separately  and  then 
transported  to  the  boiler  along  separate  lines;  as  a  result,  combus¬ 
tion  can  be  optimised  by  choosing  an  appropriate  burner  level  for 
biomass  injection. 

Case  B  can  be  further  subdivided;  in  Case  B1  the  injection  of  the 
milled  biomass  takes  place  in  an  appropriate  level  of  the  exist¬ 
ing  lignite  burner,  while  Case  B2  foresees  the  installation  of 
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dedicated  biomass  burners,  which  allow  for  further  optimisa¬ 
tion  of  the  combustion.  However,  Case  B2  requires  additional 
investment  and  increased  technical  risks;  given  the  relatively 
low  thermal  share  of  biomass  it  is  not  considered  further  in  this 
study.  Fig.  1  shows  the  different  suggested  ways  of  biomass 
feeding  and  injection  in  the  boiler  considered  in  the  study. 

3.2.  Description  of  cases  considered  for  CFD  modelling 

3.2.2.  Cases  conditions 

In  the  context  of  this  investigation,  the  fuel  quality,  the  load  of 
the  boiler  and  the  thermal  substitution  share  of  lignite  from  bio¬ 
mass  are  assumed  constant  for  all  the  cases  examined.  The  proxi¬ 
mate  and  ultimate  analysis  of  coal  and  biomass  are  given  in 
Tables  2  and  3.  The  lignite  characteristics  used  in  the  simulations 
are  based  on  the  design  lignite  data,  while  the  biomass  character¬ 
istics  are  based  on  typical  cardoon  fuel  analysis  measured  in  a 
laboratory. 

The  substitution  share  is  set  equal  to  10%  of  the  reference 
thermal  power  input;  this  value  is  chosen  by  the  plant  operator 
in  order  to  minimise  the  potential  impacts  of  the  chlorine-  and 
alkali-rich  herbaceous  biomass  combustion  in  terms  of  boiler 
slagging/fouling  and  corrosion.  In  addition,  this  respects  the  limita¬ 
tions  imposed  by  the  local  availability  of  biomass  [57],  while 
avoiding  a  very  high  specific  investment  per  unit  of  produced 
bio-energy  which  may  result  when  extensive  boiler  and/or  burner 
retrofit  is  required.  The  boiler  thermal  load  considered  in  this  work 
is  the  nominal  one. 

It  should  be  noted  that  slagging/fouling  and  corrosion  phenom¬ 
ena  is  not  examined  by  the  CFD  modelling.  Although  cardoon 
exhibits  high  chlorine,  as  well  as  alkalis,  content  that  could  cause 
such  issues  on  dedicated  combustion  units,  the  actual  chlorine 
content  in  the  fuel  mixture  for  a  10%  thermal  share  is  very  low. 
In  addition,  the  high  sulphur  content  of  lignite  results  in  a  S/Cl  mo¬ 
lar  ratio  much  higher  than  4,  which  again  indicates  non-corrosive 
behaviour  [58]. 

For  CFD  modelling  a  further  subdivision  of  the  investigated 
cases  is  based  on  the  investigation  of  (a)  the  injection  points  of  bio¬ 
mass  in  the  boiler  and  (b)  the  mean  diameter  of  biomass  particles. 

3. 2. 2. 2.  Injection  points  of  biomass.  The  insertion  position  of  bio¬ 
mass  is  associated  with  the  combustion  system  simulated.  As  pre¬ 
viously  mentioned,  in  Case  A,  the  biomass  is  fed  along  with  the 
lignite  in  the  same  already  existing  belts  and  subsequently  are 
both  pulverized  in  the  lignite  mills.  In  this  case  the  biomass  enters 
the  combustion  chamber  from  all  burners  and  all  levels,  as  lignite 
does.  The  distribution  of  cardoon  per  burner  level  is  assumed  to 
follow  that  of  lignite,  as  described  in  the  plant  operational  data. 


Fig.  1.  Cases  for  biomass  injection  in  the  boiler:  Case  A  -  “co-milling”,  Case  B1  - 
dedicated  biomass  milling/existing  lignite  burners,  Case  B2  -  dedicated  biomass 
milling/dedicated  biomass  burners. 


Table  2 

Proximate  analysis  of  lignite  and  biomass. 


Proximate 

analysis 

Lignite 

Biomass 

wt%  (as 
received) 

wt%  (after 
mills) 

wt%  (as 
received) 

wt%  (after 
mills) 

Moisture 

36.21 

12.00 

15.43 

5.00 

Ash 

27.61 

38.09 

7.30 

8.20 

Volatile 

19.01 

26.22 

68.94 

77.44 

Char 

17.17 

23.69 

8.33 

9.36 

Total 

100 

100.00 

100 

100 

NCV  (MJ/kg) 

8.06 

12.04 

16 

16.78 

Table  3 

Ultimate  analysis  of  lignite  and  biomass  (n.d.: 

not  detected). 

Ultimate  analysis  (wt%  dry,  ash  free) 

Lignite 

Biomass 

C 

61.95 

50.372 

H 

4.30 

6.742 

N 

1.65 

0.819 

O 

31.25 

41.314 

S 

0.84 

0.258 

Cl 

n.d. 

0.495 

Total 

100 

100 

Case  B1  corresponds  to  the  installation  of  a  separate  milling  system 
for  biomass  and  then  to  the  transportation  of  the  milled  biomass 
through  a  pneumatic  system  to  the  lignite  boiler.  In  this  case  the 
pulverized  biofuel  is  injected  into  the  coal  duct  downstream  of 
the  lignite  mills  and  therefore  it  enters  the  furnace  through  the 
existing  lignite  burners.  In  this  case,  it  is  possible  to  choose  the  le¬ 
vel  of  injection  into  the  coal  ducts  and  the  burner  level  through 
which  the  biomass  fuel  enters.  Specifically,  Case  Bla  includes  input 
of  entire  quantity  of  biomass  through  the  upper  main  burners  and 
Case  Bib  includes  input  of  the  entire  quantity  of  biomass  through 
the  vapour  burners.  The  input  of  biomass  from  the  level  of  the  low¬ 
er  main  burners  is  not  considered  as  a  case,  since  a  preliminary 
assessment  showed  that  input  of  biofuel  from  this  low  level  only 
will  result  in  high  losses  of  fuel  as  unburnt  material  in  the  bottom 
ash. 

3. 2. 2.2.  Mean  diameter  of  biomass  particles.  The  mean  diameter  of 
biomass  particles  is  a  key  feature  of  the  combustion  system  and 
is  strongly  related  to  the  system  specifications  for  the  biomass 
grinding  system.  In  order  to  investigate  the  effect  of  the  mean 
diameter  of  biofuel  on  its  combustion  behaviour  (temperatures, 
emissions,  char  burnout)  three  cases  have  been  considered  with 
average  equivalent  diameters  of  (i)  0.5  mm,  (ii)  1.0  mm  and  (iii) 
3.0  mm.  The  higher  biomass  particle  diameter  of  3  mm  is  not  con¬ 
sidered  for  Case  B,  since  the  installation  of  a  dedicated  biomass 
milling  system  is  expected  to  lead  to  lower  particle  sizes;  hence, 
this  case  can  be  valid  only  for  the  non-optimised  milling  of  Case 
A.  It  is  noted  that  we  have  assumed  the  average  equivalent  diam¬ 
eter  of  the  biomass  particles  in  the  simulation,  for  simplification 
reasons,  though  in  reality  biomass  particles  follow  a  size  distribu¬ 
tion.  In  addition,  this  simplification  allows  for  the  derivation  of 
clear  conclusions  about  the  effect  of  the  biomass  diameter  on  the 
operational  characteristics  of  the  plant. 

On  the  other  hand,  the  mass  weighted  distribution  of  injected 
lignite  particles  is  presented  in  Table  4.  These  values  have  been  de¬ 
rived  after  taking  samples  from  the  existing  milling  system  of  the 
plant. 

A  short  description  of  the  cases  examined  in  this  paper  is  pre¬ 
sented  in  Table  5. 
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3  A  2.  Boiler  geometry  and  boundary  conditions 

The  geometry  of  the  boiler  of  Meliti  power  plant  and  the  imple¬ 
mented  numerical  grid  are  shown  in  Figs.  1  and  2.  It  is  important  to 
note  that  the  geometry  modelled  is  in  full  accordance  with  the 
drawings  of  the  plant  without  the  usual  geometrical  simplifica¬ 
tions  in  similar  studies.  The  purpose  of  this  study  is  to  simulate 
as  accurately  as  possible  the  actual  operation  of  the  boiler,  both 
for  the  reference  and  the  co-firing  cases  (see  Fig.  3). 

The  output  of  the  furnace  and  the  flue  gas  recirculation  zones 
are  modelled  setting  boundary  conditions  of  Neumann  type  (first 
derivative  equal  to  zero).  The  output  mass  flow  of  gas  recirculation 
ducts  and  furnace  outlet  is  based  on  thermodynamic  calculations 
in  the  existing  drying  and  grinding  system,  according  to  opera¬ 
tional  data  provided  by  the  operators  of  the  plant.  Usually,  the  flue 
gas  mass  flow  rate  exiting  the  furnace  is  in-between  82.4%  and 
83.4%  of  the  total  gas  flow  inlet,  while  the  rest  is  driven  towards 
the  recirculation  ducts  in  order  to  supply  heat  for  the  fuel  drying 
in  the  mills.  The  operating  conditions  for  the  reference  case  and 
the  co-firing  ones  are  listed  in  Tables  6  and  7.  It  should  be  noted 
that  the  reference  case,  for  which  only  lignite  is  fed  to  the  boiler, 
is  based  on  the  plant  design  data  and  their  verification  by  the  plant 
acceptance  tests. 

3.2.  Description  of  cases  considered  for  economic  evaluation 
32 A.  Case  A 

As  mentioned,  Case  A  assumes  that  cardoon  is  in  the  form  of 
pellets  by  truck  from  a  long-term  storage  site  or  directly  from 
the  pellet  plant  to  the  power  plant  site.  Pellets  can  be  transported 
within  the  plant  using  the  existing  open  lignite  feeding  conveyor 
belts,  since  due  to  their  specific  weight,  pellet  losses  during  con¬ 
veying  are  not  expected  to  be  high.  In  addition,  the  existing  lignite 
mills  are  expected  to  be  able  to  handle  the  milling  of  pellets  for  the 
investigated  level  of  thermal  substitution. 

In  particular,  Case  A  requires  the  following  key  components: 

1.  The  truck  delivery  and  unloading  area  (bulk  pellets  transport). 

2.  The  intermediate  storage  system.  First  pellets  are  unloaded  in  a 
small  size  silo,  which  is  located  at  a  lower  level  for  ease  of 
unloading.  From  there,  the  pellets  are  raised  through  a  bucket 
elevator  in  a  silo  for  intermediate  storage  with  sufficient  capac¬ 
ity  to  meet  the  needs  for  48  h.  The  loading  of  silos  and  the  oper¬ 
ation  of  the  bucket  elevator  will  not  be  continuous  but 
interrupted  between  each  new  delivery  (1  truck  per  hour 
estimated). 

3.  The  feeding  system  from  the  silo  to  the  main  belt-mixing  with 
lignite.  From  the  intermediate  storage  silo  the  pellets  are  driven 
through  a  screw  feeder  at  a  weighting  belt  and  from  there 
through  a  conveyor  belt  at  a  distance  of  100  m  from  the  silo 
to  the  main  conveyor  belt  where  it  is  mixed  with  lignite.  The 
mixture  is  led  to  the  coal  crushers  and  then  the  coal  silo.  From 
there  it  is  fed  to  the  mills  and  then  to  existing  lignite  burners. 

The  following  data  are  considered: 

•  Daily  consumption  of  biomass:  500  t/day. 

•  Bulk  density  of  pellets:  600  kg/m3. 

•  Required  silo  storage  capacity:  2  *  500/0.6  =  1666.7  m3. 


Table  4 

Mass  weighted  diameter  distribution  for  lignite  particles. 


Diameter  (pm) 

Mass  percentage  w.w% 

300 

39.8 

170 

14.01 

90 

12.08 

50 

34.11 

Table  5 

Short  description  of  the  cases  examined. 


Biomass  particle  size 
distribution 

Injection  methods/cases 

A 

B1  (a)  -  upper  main 
burners 

B1  (b)  -  vapour 
burners 

i  (0.5  mm) 

X 

X 

X 

ii  (1.0  mm) 

X 

X 

X 

iii  (3.0  mm) 

X 

•  Loading  of  silos  and  bucket  elevator  operation:  interrupted, 
after  each  new  delivery  of  truck  (1  truck/h  estimated). 

•  Silo  unloading  and  operation  of  feeding  system  and  conveyor 
belt:  interrupted  according  to  lignite  loading. 

•  Lignite  loading  duration:  10  h  per  day. 

•  Mass  flow  of  biomass  in  the  feeding  system  from  the  feeding 
silo:  500 1/10  h  =  50  t/h  or  50/0,  6  =  83.3  m3/h. 

3.2.2.  Case  B 

Case  B  assumes  that  the  biomass  fuel  is  delivered  as  bales  (50%) 
or  as  shredded  material  (50%).  The  existing  handling  system  for  lig¬ 
nite  cannot  be  used  and  there  is  a  need  for  dedicated  milling  and 
conveying  systems. 

In  this  case,  an  unloading  and  delivery  area  for  trucks  transport¬ 
ing  the  biomass  in  bales  or  in  shredded  bulk  form  is  also  provided. 
The  case  of  receiving  50%  of  the  incoming  quantity  of  biomass  into 
bales  and  50%  in  bulk  form  is  considered.  The  biomass  bales  are 
now  stored  in  a  covered  space,  protected  from  weather  conditions, 
in  which  there  is  the  capability  of  loading  them  with  a  loader,  and 
an  automated  bag  opening  system  for  the  bales.  Similarly,  in  the 
case  of  receiving  shredded  biomass  in  bulk  form,  it  will  be  stored 
in  piles  inside  the  storage  building  or  in  covered  silos  with  auto¬ 
matic  walking  floor  feeding  system. 

In  particular,  this  case  includes  the  following  components: 

1.  Truck  delivery  and  unloading  area. 

2.  The  intermediate  storage  system  (area  near  the  station 
entrance).  The  biomass  bales  are  stored  in  a  covered  space,  pro¬ 
tected  from  weather  conditions,  in  which  there  is  the  capability 
of  loading  them  with  a  loader,  and  an  automated  bag  opening 
system  for  the  bales.  Similarly,  in  the  case  of  receiving  shredded 
biomass  in  bulk  form,  it  will  be  stored  in  piles  inside  the  storage 
building  or  in  covered  silos  with  automatic  walking  floor  feed¬ 
ing  system.  The  storage  system  should  be  large  enough  to 
accommodate  enough  biomass  for  2  days  of  operation  (1000 1). 

3.  Pretreatment,  feeding  and  grinding  of  biomass.  From  the  inter¬ 
mediate  storage,  biomass  is  transferred  to  a  separate  biomass 
grinding  mill  where  it  is  pulverized  and  powdered  biomass  is 
stored  temporarily  and  driven  pneumatically  up  to  the  combus¬ 
tion  system. 

4.  Buffer  storage,  pneumatic  feeding  to  the  boiler.  From  this  inter¬ 
mediate  storage,  biomass  is  transferred  to  a  separate  biomass 
grinding  mill  where  it  is  pulverized  and  powdered  biomass  is 
stored  temporarily  and  driven  pneumatically  up  to  the  combus¬ 
tion  system.  Both  the  buffer  silo  and  the  pneumatic  feeding  sys¬ 
tems  should  follow  the  APEX  standards. 

The  following  data  are  taken  into  account: 

•  Daily  consumption  of  biomass:  500  t/days. 

•  Mass  flow  of  biomass  from  storage  to  mill  grinding:  continuous 
500  t/d/24  h  =  20.8  t/h. 

•  Biomass  mill  operation:  continuous. 

•  Intermediate  storage  of  powder:  buffer  for  a  short  time  1-2  h. 

•  Pneumatic  biomass  powder  feeding  up  to  the  combustion  sys¬ 
tem:  continuous  20.8  t/h,  at  a  distance  of  350  m. 
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Fig.  2.  (a)  Geometry  of  Meliti  Boiler,  (b)  top  view  of  Meliti  Boiler  and  (c)  numerical  grid. 


re-bumers 


secondary 


VaP°ur  Lower  main  S  Upper  main 

burner  burner  burner 


Fig.  3.  (a)  Over-fire  air  (OFA)  nozzle  and  recirculation  ducts  geometry  and  (b)  primary  and  secondary  airstream  in  the  main  burners,  oil  burners  and  vapour  burners 
geometry. 


4.  Results  and  discussion 

4.1.  Results  of  technical  evaluation  (CFD  analysis) 

4.1.1.  Reference  conditions 

For  the  reference  case,  lignite  is  the  only  solid  fuel  entering  the 
boiler;  it  is  injected  into  six  (6)  out  of  the  total  eight  (8)  burners. 


The  main  furnace  parameters  for  the  reference  case  are  presented 
in  Table  8  and  compared  to  design  values  or  data  from  the  plant 
acceptance  tests.  More  specifically,  the  furnace  outlet  temperature 
is  comparable  to  the  design  value  of  992  °C,  while  the  total  heat 
flux  through  the  membrane  walls  was  found  to  be  413  MWth.  This 
is  in  good  agreement  with  a  heat  transfer  of  396.9  MWth  estimated 
by  taking  into  account  the  following  design  data  and  estimates: 
steam  mass  flow  (282.4  kg/s),  inlet  and  outlet  steam  data 
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Table  6 

Boiler  operating  conditions  for  the  reference  case. 


Stream 

Mass 

flow 

(kg/s) 

Temperature 

(°C) 

Gas  composition  (wt%) 

O2  H2O  CO2 

n2 

Primary  air 

283.54 

220 

13.55 

16.03 

6.83 

63.59 

(transport  air) 
Secondary  air 

164.10 

305 

23.21 

76.79 

(from  boilers) 
Secondary  air  (by 

66.37 

305 

23.21 

. 

. 

76.79 

stage 

combustion 
system) 
Secondary  air 

11.52 

151 

23.21 

76.79 

(from  ash 
hopper) 

Lignite 

117.0 

100 

Biomass 

— 

— 

— 

— 

— 

— 

Table  7 

Boiler  operating  conditions  for  co- 

-firing  with  10%  biomass  thermal  share. 

Stream 

Mass 

flow 

(kg/s) 

Temperature 

(°C) 

Gas  composition  (wt%) 

O2  H2O  CO2 

n2 

Primary  air 

272.92 

220 

13.91  15.26 

6.63 

64.20 

(transport  air) 
Secondary  air 

162.83 

305 

23.21  - 

76.79 

(from  boilers) 
Secondary  air  (by 

65.86 

305 

23.21  - 

. 

76.79 

stage 

combustion 
system) 
Secondary  air 

11.43 

151 

23.21  - 

76.79 

(from  ash 
hopper) 

Lignite 

105.3 

100 

Biomass 

6.42 

100 

-  - 

- 

- 

(316  °C/ 268.7  bar  and  418  °C/ 258.1  bar  respectively)  and  3%  radia¬ 
tion  losses  for  the  total  steam  production  of  the  boiler  (870  MWth). 
The  02  and  NO*  concentration  values  for  the  reference  case  are  also 
in  line  with  the  plant  acceptance  tests  and  the  typical  values  of 
operation  under  full-load  conditions  (02:  2.9  wt%  compared  to 
2.79  wt%,  NO*  ~  254  mg/N  m3  compared  to  200  mg/N  m3,  both  at 
6  vol.%  dry  02).  From  this  comparison,  it  is  evident  that  the  numer¬ 
ical  model  predicts  quite  well  the  main  operational  characteristics 
of  the  power  plant,  thus  allowing  us  to  proceed  with  the  examina¬ 
tion  of  the  co-firing  concepts. 

4.1.2.  Co -firing  simulations  -  investigation  of  the  biomass  particle  size 
for  Case  A 

The  basic  operating  parameters,  as  obtained  numerically  for 
cases  A  (i),  A  (ii)  and  A  (iii)  are  also  presented  and  compared  to 
the  reference  case  in  Table  8.  Based  on  the  results,  no  significant 
operational  changes  are  expected  to  occur  in  the  boiler  character¬ 
istics,  except  for  an  increase  of  around  4%  in  the  thermal  flux  to¬ 
wards  the  membrane  walls.  The  concentration  of  NOx  varies  by 
up  to  8%  with  a  decreasing  trend,  showing  that  the  recovery  of  lig¬ 
nite  with  biomass  has  a  positive  impact  on  NO*  emissions.  Table  9 
shows  the  efficiency  of  the  furnace  that  is  calculated  based  on  the 
wet  and  solid  ash  losses. 

Fig.  4,  depicts  the  distribution  of  the  mean  temperature  along 
the  height  of  the  boiler  being  developed  for  the  reference  and  cases 
A  (i),  A  (ii)  and  A  (iii).  As  can  be  seen,  the  mean  temperature  is  prac¬ 
tically  the  same  for  all  examined  cases  for  heights  higher  than 
20  m.  At  lower  heights,  higher  temperatures  are  observed 


compared  to  the  reference  case  as  the  biomass  particle  size  in¬ 
creases;  the  mean  temperature  of  Case  A  (iii)  at  a  height  of  12  m 
is  500  °C  higher  compared  to  the  reference  case.  This  is  due  to 
the  shift  of  the  volatile  release  and  combustion  section  of  the  fur¬ 
nace  to  lower  heights  as  the  biomass  particle  size  increases.  These 
local  hot  spots  could  cause  problems  in  the  boiler  operation  if  the 
substitution  share  is  much  higher  than  10%  and  should  be  studied 
in  depth  in  the  future. 

The  total  lignite  and  biomass  char  burnout  are  respectively  pre¬ 
sented  in  Figs.  5a  and  5b,  both  for  the  whole  volume  of  the  boiler 
and  for  the  particles  exiting  the  furnace  top  and  hopper.  Further¬ 
more,  Figs.  6a  and  6b  present  the  mass  flow  percentage  of  lignite 
and  biomass  particles  exiting  the  boiler  outlet  and  hopper  for  each 
burner  level.  The  rest  of  the  fuel  particles  are  exiting  the  boiler 
through  the  recirculation  ducts. 

From  these  graphs  the  following  conclusions  can  be  extracted: 

•  The  lignite  and  biomass  entering  from  vapour  burners  (re-burn- 
ers)  exit  through  the  furnace  outlet  with  a  percentage  of 
approximately  87%,  while  the  rest  passes  through  the  recircula¬ 
tion  ducts  regardless  of  the  particle  size  of  biomass. 

•  The  percentage  of  the  lignite  particles  from  the  lower  burners 
that  exit  the  furnace  through  the  hopper  increases  to  28.4% 
for  Case  A  (i)  and  29.08%  for  Case  A  (ii)  compared  to  26.3%  for 
the  reference  case;  a  decrease  to  24.47%  is  observed  for  Case 
A  (iii).  For  cases  A  (i)  and  A  (ii)  it  appears  that  the  lower  flue 
gas  velocities  in  the  hopper  allow  for  increased  lignite  particle 
losses;  for  Case  A  (iii)  the  high  temperature  increase  in  the  hop¬ 
per  area  compensates  for  the  reduced  flue  gas  flow  and  results 
in  higher  velocities  and,  therefore,  to  a  reduction  of  lignite  par¬ 
ticles  that  fall  on  the  ash  hopper. 

•  The  corresponding  percentage  for  biomass  leaving  the  ash  hop¬ 
per  increases  significantly  for  Case  A  (iii)  (Reference  case:  0%, 
Case  A  (i):  0.21%,  Case  A  (ii):  0.74%,  Case  A  (iii):  31.25%),  sug¬ 
gesting  that  the  diameter  of  3  mm  presents  significant  ash  hop¬ 
per  losses  and  therefore  should  be  avoided  as  an  operating  case. 

•  Note  that  the  biomass  char  burnout  is  significantly  reduced  as 
the  diameter  increases,  particularly  for  biomass  particles  enter¬ 
ing  from  the  upper  and  lower  main  burners. 

From  these  results  one  can  conclude  that  Case  A  (iii)  should  be 
avoided  since  it  results  in  significant  mass  and  energy  losses  of  so¬ 
lid  particles  in  the  bottom  ash,  while  these  losses  are  significantly 
greater  for  the  lower  burners  in  comparison  with  the  upper. 


Table  8 

Comparison  of  main  furnace  parameters  between  design  data,  reference  and  co-firing 
cases. 


Case 

Temperature 

(°C) 

O2 

(wt%) 

co2 

(wt%) 

Total  heat 

flow 

through 

membrane 

walls 

(MWth) 

NOx 

mg/ 

Nm3 
(6  vol.% 
dry) 

Design  data/ 

992 

2.79 

20.0 

396.9 

<=200 

acceptance 

tests 

Reference 

983 

2.90 

19.48 

412.97 

253.8 

A  (i)  -  0.5  mm 

993 

3.32 

19.55 

434.60 

234.8 

A  (ii)  -  1.0  mm 

987 

3.46 

19.40 

427.90 

242.9 

A  (iii)  -  3.0  mm 

985 

3.52 

19.40 

428.55 

233.86 

Bla  (i) 

995 

3.08 

19.82 

443.19 

257.96 

Bib  (i) 

1005 

3.08 

19.83 

441.18 

212.76 

Bla  (ii) 

991 

3.14 

19.74 

443.78 

246.2 

Bib  (ii) 

1000 

3.14 

19.74 

434.97 

231.2 
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Table  9 

Furnace  performance  based  on  combustion  efficiencies  for  ash  losses. 


Combustion 
efficiency 
based  on  fly 
ash  losses  (%) 

Combustion 
efficiency 
based  on  wet 
ash  losses  (%) 

Weighted  combustion 
efficiency  based 
on  wet  and 
fly  ash  losses  (%) 

Reference 

99.94 

84.24 

97.32 

A  (i)  -  0.5  mm 

99.90 

89.14 

98.13 

A  (ii)  -  1.0  mm 

99.86 

85.57 

97.41 

A  (iii)  -  3.0  mm 

99.13 

84.56 

96.96 

Bla  (i) 

99.96 

93.34 

99.07 

Bib  (i) 

98.49 

94.50 

98.11 

Bla  (ii) 

99.95 

93.14 

99.00 

Bib  (ii) 

98.03 

91.27 

97.11 

Fig.  7  shows  the  contours  of  temperature  in  the  vertical  plane 
passing  crossing  the  centre  of  the  boiler  (X  =  0.0),  and  the  total  heat 
flux  on  membrane  walls.  Fig.  8  presents  the  ratio  of  the  radiative 
heat  flux  to  the  total  heat  flux;  as  can  be  seen,  the  heat  leaving 
the  furnace  through  the  membrane  walls  is  mostly  due  to  radia¬ 
tion,  as  is  expected  of  such  boiler  types.  Fig.  9  shows  the  tempera¬ 
ture  and  velocity  distributions  and  concentrations  of  NO*  at  the 
burner  levels. 

Therefore  a  10%  thermal  substitution  share  does  not  signifi¬ 
cantly  affect  the  proper  operation  of  the  boiler  in  terms  of  the 
developing  temperature  fields  and  NOx  emissions,  except  for  the 
Case  A  (iii)  (3  mm),  where  biomass  losses  from  ash  hopper  are 
quite  high. 

4.1.3.  Co -firing  simulations  -  investigation  of  the  biomass  particle  size 
for  injection  at  different  burner  levels  (Cases  Bla  and  Bib) 

4.13.1.  Case  (i)  -  0.5  mm:  Cases  Bla  (i)  and  Bib  (i).  Table  8  presents 
the  basic  operating  parameter  values  for  these  cases  compared  to 
the  reference  one.  From  the  table  it  is  shown  that  the  boiler  oper¬ 
ation  is  not  significantly  affected  in  any  case  with  10%  substitution 
of  thermal  power  from  biomass,  except  for  the  fact  that  there  is  an 
increase  in  the  total  membrane  wall  heat  flux  of  about  7%  and  a 


significant  NO*  reduction  by  around  16.2%  in  the  case  of  injecting 
biomass  from  the  vapour  burners. 

In  these  cases,  the  mixed  combustion  improves  the  efficiency  of 
the  furnace.  Fig.  10  depicts  the  mean  temperature  spatial  distribu¬ 
tion  in  the  boiler  for  reference  Case  A  (i),  Bla  (i)  and  Bib  (i). 

For  these  cases  no  significant  changes  in  the  boiler  operation 
are  expected  when  biomass  enters  the  furnace  through  individ¬ 
ual  fuel  inlets  (upper  main  and  vapour  burners).  As  the  height 
of  biomass  injection  increases,  the  position  of  the  hot  spots  are 
developing  in  higher  heights  (the  difference  in  temperature  at 
12  m  for  cases  A  (i)  and  Bib  (i)  is  about  550  °C).  Thus,  the 
biomass  injection  from  one  burner  level  is  not  optimal  for  main¬ 
taining  the  flame  temperature  uniform  throughout  the  entire 
furnace  height  range,  especially  if  the  thermal  share  substitution 
is  much  higher  than  10%.  In  Figs.  11a  and  lib,  the  total  lignite 
and  biomass  char  burnouts  are  presented.  Furthermore, 
Figs.  12a  and  12b  present  the  mass  flow  percentage  of  lignite 
and  biomass  solid  particles  exiting  the  boiler  outlet  and  hopper 
per  burner  level. 

From  the  afore-cited  graphs  the  following  conclusions  can  be 
extracted: 

•  For  cases  where  biomass  enters  the  furnace  in  specific  burners 
(either  the  upper  main  ones  or  the  vapour  burners)  the  lignite 
char  burnout  is  improving  compared  to  the  case  of  biomass 
entering  from  all  the  three  burner  levels. 

•  In  each  of  the  cases  studied,  the  amount  of  lignite  leaving  the 
ash  hopper  is  decreasing,  when  it  enters  from  the  lower  main 
burners  (Case  A  (i):  28.41%,  Case  Bla  (i):  23.82%,  Case  Bib  (i): 
21.47%)  and  upper  main  burners  (Case  A  (i):  4.19%,  Case  Bla 
(i):  3.16%,  Case  Bib  (i):  1.33%),  improving  the  combustion  qual¬ 
ity  of  the  boiler.  This  suggests  reduced  possibility  for  unburnt 
lignite  losses  in  the  hopper. 

•  In  the  case  when  biomass  enters  through  the  upper  main  burn¬ 
ers  only  (Bla),  there  is  a  little  increase  in  the  amount  of  biomass 
leaving  the  outlet,  while  if  the  biomass  comes  from  vapour 
burners  only  (Bib)  the  amount  is  reduced.  For  both  cases  the 
ash  hopper  losses  are  minimal  for  biomass. 


Fig.  4.  Distribution  of  mean  temperature  along  the  boiler  height  for  the  reference  case,  A  (i),  A  (ii)  and  A  (iii). 
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□  Lignite -Total  □  Lignite  -  Outlet  □  Lignite  -  Hopper 
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Fig.  5a.  Lignite  char  burnout  per  burner  level,  overall,  at  output  and  at  ash  hopper  of  the  boiler. 
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Fig.  5b.  Biomass  char  burnout  per  burner  level,  overall,  at  output  and  at  ash  hopper  of  the  boiler. 


To  conclude,  the  combustion  behaviour  of  biomass  is  not 
changing  significantly,  while  the  behaviour  of  lignite  both  in 


terms  of  burnout  and  ash  hopper  losses  is  improving,  making 
the  Bib  case  the  best  to  choose  among  all  the  examined.  The 
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□  Lignite -Outlet  □  Lignite  -  Hopper 


□  Biomass  -  Outlet  ■  Biomass  -  Hopper 


selection  of  Case  Bib  as  the  optimal  one  is  also  supported  by  the 
significant  reduction  in  the  NOx  concentration.  Despite  the 


observed  temperature  increases  and  the  expectations  for  in¬ 
creased  thermal  NOx  formation,  the  main  reason  for  the  decrease 
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Fig.  7.  (a)  Contours  of  temperature  on  the  vertical  plane  and  (b)  heat  flux  on  walls,  for  cases  A. 


of  the  total  NO*  emissions  during  biomass  co-firing  is  the  high  de¬ 
crease  of  fuel  NO*  due  to  the  different  intermediate  species  paths 
for  the  two  fuels.  As  discussed  in  detail  in  Section  2.1.3,  NH3  and 
HCH  are  the  main  intermediate  species  during  biomass  and  lig¬ 
nite  N-conversion  respectively;  the  former  leads  to  lower  fuel 
NO*  over  the  latter. 


4.1.32.  Case  (ii)  -  1.0  mm:  Cases  Bla  (ii)  and  Bib  (ii).  The  mean 
diameter  in  this  investigation  is  1.0  mm.  Tables  8  and  9  present 
the  basic  operating  conditions  and  furnace  performance  for  cases 
Bla  (ii)  and  Bib  (ii)  compared  to  the  reference  one.  Again  the  boi¬ 
ler  operation  is  not  significantly  affected,  except  for  the  fact  that 
there  is  an  average  increase  in  total  heat  flux  through  membrane 
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Fig.  7.  ( continued ) 


walls  by  approximately  6.5%  and  a  significant  reduction  of  up  to 
8.7%  in  the  NO*  emissions  when  increasing  the  biomass  injection 
height.  The  increased  heat  flux  towards  the  membrane  walls  is 
justified  by  the  lower  flue  gas  mass  flow  rate  for  the  co-firing 
cases  compared  to  the  reference  case,  coupled  with  the  fact  that 
the  flue  gas  temperature  at  the  outlet  of  the  furnace  is  not 
increasing  significantly.  Therefore,  the  enthalpy  of  the  flue  gas  at 
the  furnace  exit  is  lower  and  since  the  thermal  input  remains 
the  same,  the  energy  balance  requires  that  the  heat  flux  towards 
the  membrane  walls  is  higher  in  the  co-firing  cases. 


In  these  cases  mixed  combustion  also  improves  the  efficiency 
of  the  furnace.  Fig.  13  shows  the  distribution  in  height  of  the 
average  temperature  in  the  boiler  being  developed  for  the 
reference  operation,  Case  A  (ii),  Bla  (ii)  and  Bib  (ii).  The  afore¬ 
cited  graph  reveals  that  the  induced  temperature  field  follows 
the  same  behaviour  with  the  corresponding  of  cases  (i).  The  ten¬ 
dency  of  the  phenomenon  to  present  a  delayed  development  of 
hot  temperature  spots  is  not  significantly  affected  and 
differentiated  between  case  (i)  (0.5  mm)  and  case  (ii) 

(1.0  mm). 
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Fig.  8.  Radiative  heat  flux  to  total  heat  flux  on  membrane  walls  for  cases  A. 


Figs.  14a  and  14b  present  the  total  lignite  and  biomass  char 
burnouts.  Furthermore,  Figs.  15a  and  15b  present  the  mass  flow 
percentage  of  lignite  and  biomass  solid  particles  exiting  the  boiler 
outlet  and  hopper  per  burner  level. 

From  these  graphs  the  following  important  conclusions  can  be 
extracted: 


•  For  cases  where  the  biomass  enters  only  from  the  upper  main 
burners,  an  improved  lignite  char  burnout  is  clear  compared 
to  the  case  of  biomass  entering  from  all  the  three  burner  levels. 
When  the  biomass  enters  through  the  vapour  burners,  lignite 
burnout  slightly  decreases. 
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Fig.  9.  (a)  Induced  velocity  field  (m/s)  and  (b)  NOx  concentration  spatial  distribution  at  burner  levels  for  cases  A. 


•  The  amount  of  lignite  leaving  the  ash  hopper  is  reduced  when 
entering  from  the  lower  main  burners  (Case  A  (ii):  29.08%,  Case 
Bla  (ii):  24.60%,  Case  Bib  (ii):  24.20%)  and  upper  main  burners 
(Case  A  (ii):  4.1 7%,  Case  B1  a  (ii):  3.74%,  Case  Bib  (ii):  2.63%).  This 
suggests  reduced  possibility  for  unburnt  lignite  losses  in  the 
hopper. 

•  In  Case  Bla  (ii)  (entrance  from  upper  main  burners),  there  is  lit¬ 
tle  increase  in  the  amount  of  biomass  leaving  the  outlet  but  sig¬ 
nificant  reduction  in  the  burnout  of  biomass  exiting  the  hopper 
(Case.  A  (ii):  100%,  Case  Bla  (ii):  93.43%).  In  Case  Bib  (entrance 
from  vapour  burners),  the  amount  of  biomass  leaving  the  outlet 


decreases  (Case  A  (ii):  86.88%,  Case  Bib  (ii):  84.11%).  For  both 
cases  the  ash  hopper  losses  are  small  for  biomass,  just  as  in 
the  previous  operating  cases  (i). 

To  conclude,  the  combustion  behaviour  of  biomass  in  these 
cases  does  not  change  significantly,  while  the  behaviour  of  lignite 
both  in  terms  of  burnout  and  ash  hopper  losses  are  highly  im¬ 
proved.  For  Case  Bla  (ii)  however,  where  biomass  enters  through 
the  upper  main  burners,  ash  hopper  losses  are  increased  by  0.07% 
along  with  a  significant  decrease  in  the  biomass  burnout.  Therefore 
Case  Bib  (ii)  is  considered  preferable  compared  to  Bla  (ii). 
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Fig.  9.  ( continued ) 


All  the  aforementioned  predictions  indicate  that  biomass 
should  be  entering  the  boiler  with  a  mean  diameter 
between  0.5  and  1  mm,  so  that  combustion  without  significant 
losses  can  be  maintained.  Nevertheless,  the  most  optimal  selec¬ 
tion  of  operation  is  for  biomass  to  be  injected  from  the  vapour 
burners. 


4.2.  Results  of  economic  evaluation 

4.2. 1 .  Estimation  of  expected  costs  and  revenues  from  the  application 
of  co-firing  cases 

The  difference  in  Net  Present  Value  of  each  investment  in  rela¬ 
tion  to  the  case  of  non-mixed  combustion  is  affected  by: 


•  Additional  capital  and  operating  expenses  due  to  the  operation 
of  co-firing  systems.  These  costs  include:  (a)  Investment  cost, 

(b)  The  difference  in  fuel  costs  due  to  the  change  of  fuel,  and 

(c)  the  additional  costs  of  operation  and  maintenance,  person¬ 
nel,  etc. 

•  The  additional  revenue/savings  generated  by  biomass  co-firing. 
It  includes:  (a)  the  increased  revenues  from  electricity  genera¬ 
tion  due  to  the  biomass  feed-in  tariff  and  (b)  the  savings  from 
C02  emission  savings. 


4.2.1.  t .  Total  cost  of  investment.  The  Total  Plant  Cost  (TPC)  includes 
design  costs,  procurement  costs  and  construction  costs.  For 
co-firing,  a  wide  range  of  100-600  USD/kWe  is  reported  [59], 
depending  on  the  fuel  and  technical  solutions  adopted.  Projects 
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Fig.  10.  Distribution  of  average  temperature  of  the  boiler  in  height  for  the  reference  case,  A  (i),  Bla  (i)  and  Bib  (i). 


□  Lignite -Total  ■  Lignite  -  Outlet  □  Lignite  -  Hopper 


Fig.  11a.  Lignite  burnout  per  burner  level,  overall,  at  furnace  outlet  and  at  ash  hopper. 


that  implement  modest  technological  solutions  and  minimal  mod¬ 
ifications  in  the  boiler  section  require  much  lower  specific  costs 
compared  to  ones  where  extensive  boiler  retrofitting  is  required. 
In  any  case,  these  values  are  much  lower  compared  with  the  corre¬ 
sponding  specific  costs  of  new  facilities  of  100%  dedicated  biomass 


combustion,  which  may  reach  2500-3000  €/kWe.  This  significant 
difference  in  cost  is  due  to  existing  equipment  used  in  the  case 
of  co-firing. 

In  this  work,  the  required  TPCs  for  the  implementation  of 
Cases  A  and  B  is  based  on  economic  offers  received  by  the  plant 
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Fig.  lib.  Biomass  char  burnout  per  burner  level,  overall,  at  furnace  outlet  and  at  ash  hopper. 


□  Lignite  -  Outlet  ■  Lignite  -  Hopper 


operator  for  the  technical  specifications  described  in  Section  3.2. 
The  economic  evaluation  is  based  on  the  Total  Cost  of  Investment 
(TCI),  also  known  as  the  Capital  Expenditure  (CAPEX),  which  is 
calculated  by  increasing  the  TPC  by  5%  for  development  costs 


and  by  a  further  10%  as  additional  unplanned  costs.  The  TCI  is 
given  in  Table  12  and  amounts  to  15.6€/kWe  for  Case  A  and 
41.7  €/l<We  for  Case  B.  As  explained,  these  low  values  are  justified 
by  the  lack  of  any  boiler  retrofitting.  It  should  be  noted  that  the 
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Fig.  12b.  Mass  percentage  of  biomass  particles  leaving  the  outlet  and  the  ash  hopper  of  the  boiler  per  burner  level. 


Fig.  13.  Distribution  of  average  temperature  in  height  of  the  boiler  for  the  reference  case,  A  (ii),  Bla  (ii)  and  Bib  (ii). 


70%  of  the  TCI  is  considered  to  be  covered  by  own  funds  of  the 
operator,  while  the  remaining  30%  comes  from  long-term  loans, 
with  an  interest  rate  of  7%. 

4.2. 1.2.  Operation  and  maintenance  cost.  The  cost  of  Operation  and 
Maintenance  (O&M)  is  divided  into  fixed  and  variable.  The  fixed 


includes  all  fixed  costs  that  do  not  depend  on  annual  production, 
such  as  personnel  costs,  security  and  cost  of  planned  mainte¬ 
nance  work.  The  variable  costs  represent  the  costs  arising  from 
the  operation  of  the  unit  mainly  including  electricity  consump¬ 
tion  and  consumables.  In  this  work,  the  O&M  cost  is  considered 
as  35%  of  the  total  investment  cost,  based  on  initial  cost 
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□  Lignite -Total  ■  Lignite  -  Outlet  □  Lignite  -  Hopper 


Fig.  14a.  Lignite  burnout  per  burner  level,  overall,  at  furnace  outlet  and  at  ash  hopper. 
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Fig.  14b.  Biomass  burnout  per  burner  level,  overall,  at  furnace  outlet  and  at  ash  hopper. 


estimates  of  the  electricity  consumption,  man-hours  and  spare- 
part  costs  for  each  case. 


42 A 3.  Fuel  cost.  Co-firing  requires  significant  fuel  costs  for  the 
purchase  of  biomass  fuel.  Currently,  there  is  no  established 
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□  Lignite  -  Outlet  a  Lignite  -  Hopper 
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Fig.  15b.  Mass  percentage  of  biomass  particles  leaving  the  outlet  and  the  ash  hopper  of  the  boiler  per  burner  level. 


biomass  market  for  co-firing  plants  in  Greece,  hence  there  is  great 
uncertainty  in  the  actual  fuel  cost;  estimations  range  from  50  €/t  to 


150€/t  depending  among  others  on  the  form  of  the  delivered 
material  (pelletized  biomass  costs  more  than  biomass  bales).  A 
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Table  10 

Summary  of  cost  categories  for  Cases  A  and  B. 


Units 

Case  A 

Case  B 

Total  investment  cost 

€ 

515,902 

1,376,294 

Additional  fuel  cost  (biomass  cost  minus 
lignite  saving) 

€/a 

17,398,575 

6,894,154 

O&M  costs 

€/a 

202,161 

620,211 

Table  11 

Additional  revenues/savings  generated  by  co-firing. 

Additional  revenue  from  biomass  feed-in  tariff  (€) 

19,767,600 

Savings  from  C02  emissions  reduction  (€) 

9,161,155 

Table  12 

Parameters  used  in  the  economic  analysis. 

Parameters 

Units 

Case  A 

Case  B 

Thermal  substitution 

% 

10 

10 

Expected  life  time 

Years 

10 

10 

Development  and  construction  period 

Years 

1 

1 

Availability 

h/year 

7200 

7200 

Biomass  cost 

€/t 

150 

80 

Lignite  cost 

€/t 

19 

19 

Average  electricity  sale  price 

€/ 

MW  h 

55 

55 

Feed-in  tariff  for  biomass  use 

€/ 

MW  h 

150 

150 

C02  emission  allowance  purchase  Cost 

e/t  co2 

15 

15 

Money  market  rate 

O/ 

Zo 

7 

7 

Percentage  of  investment  covered  by  own  funds 

O/ 

Zo 

70 

70 

Percentage  of  investment  covered  by  loans 

O/ 

Zo 

30 

30 

Borrowing  interest  rate 

O/ 

Zo 

7 

7 

Payback  period  of  loan 

Years 

10 

10 

Inflation 

o/ 

Zo 

5 

5 

Table  13 

Financial  evaluation  criteria  for  Cases  A  and  B. 


Units 

Case  A 

Case  B 

Net  Present  Value  (NPV) 

€ 

23,735,926 

89,873,631 

Internal  Rate  of  Return  (IRR) 

% 

601,071% 

998,990% 

Payback  period 

Years 

1 

1 

very  high  fuel  cost  can  threaten  the  economic  viability  of  the 
investment,  if  the  expenditure  is  much  higher  than  the  corre¬ 
sponding  revenue. 

The  price  of  lignite  is  also  important  in  the  calculation  of  the 
fuel  cost,  since  the  cost  savings  achieved  by  the  biomass  substitu¬ 
tion  must  be  taken  into  account. 

For  biomass,  an  estimated  price  of  150  €/t  is  used  as  a  reference 
for  Case  A  (pellets  obtained  from  the  market),  while  for  Case  B  (bio¬ 
mass  bales  and  shredded  material)  the  cost  reference  is  taken  as 
80  €/t.  The  cost  of  lignite  is  2.26  €/GJ  or  19  €/t. 

The  cost  breakdown  for  of  the  investigated  cases  is  presented  at 
Table  10. 

4.2  A  A.  Revenues  and  savings.  The  expected  additional  revenues 
from  the  application  of  substitution  cases  and  biomass  co-firing 
are  derived  from  the  feed-in  tariff  of  electricity  from  biomass  as 
set  by  Law  GG  3851/2010  of  the  Greek  legislation.  For  a  biomass 
capacity  higher  than  5  MWe,  as  is  the  case  for  the  investigated 
cases,  the  feed-in  tariff  it  is  set  at  150  €/MWhe. 


Thus,  the  additional  revenues  from  the  sale  of  electricity  is 
resulting  from  the  difference  in  selling  price  of  electricity  from  bio¬ 
mass  in  relation  to  System  Marginal  Price  (SMP).  The  SMP  is  actu¬ 
ally  determined  by  the  market  pool  and  has  significant  variations; 
for  the  purposes  of  this  study,  an  indicative  constant  average  value 
of  55  €/MWhe  is  considered. 

The  economic  calculations  should  also  consider  the  savings  for 
the  purchase  of  C02  emissions  allowances.  This  is  a  variable  cost 
which  can  be  reduced  through  co-firing  since  biomass  is  consid¬ 
ered  carbon-neutral.  After  the  expiry  of  the  period  of  the  Kyoto 
Protocol  in  2012,  the  entire  amount  of  C02  emissions  from  lignite 
combustion  will  be  included  in  the  emissions  trading  scheme 
(ETS)  and  will  have  to  be  purchased  in  the  market.  It  is  therefore 
expected  that  the  operating  costs  of  outdated  lignite  stations  with 
low  efficiencies  will  be  significantly  affected.  In  this  context,  co-fir¬ 
ing  is  a  realistic  and  reliable  option  that  can  bring  additional  reve¬ 
nue  or  potential  savings  of  future  costs  for  existing  lignite-fired 
plants.  For  the  emissions  trading  market,  a  typical  price  of  15€/t 
C02  is  considered,  which  corresponds  approximately  to  the  aver¬ 
age  of  recent  years.  Studies  suggest  an  expected  future  increase 
of  emissions  allowance  prices  at  around  22  €/t  C02  in  the  coming 
years.  Nevertheless,  the  calculations  employed  the  modest  price 
of  15  €/t  C02  and  the  effect  of  variations  is  studied  through  para¬ 
metric  risk  analysis. 

The  additional  revenues  and  savings,  as  calculated,  are  given  in 
the  Table  11  (common  to  both  cases). 

Apart  from  the  revenues  and  expenses  which  were  mentioned 
above,  from  which  the  annual  cash  flow  of  the  investment  is  calcu¬ 
lated,  in  order  to  sum  up  the  Net  Cash  Flow  for  the  calculation  of 
the  Net  Present  Value,  amortization  instalments  required  to  repay 
the  loan  must  be  also  taken  into  account.  An  amount  attributable 
to  taxes  must  also  be  deducted,  assuming  a  tax  rate  equal  to  25% 
of  income. 

The  major  economic  parameters  used  in  the  cases  considered 
for  the  economic  analysis  are  shown  in  Table  12. 

Based  on  the  aforementioned  parameters  and  values,  the  corre¬ 
sponding  values  for  NPV,  IRR  and  payback  period  for  Cases  A  and  B 
are  presented  in  Table  13.  While  both  Cases  are  profitable,  Case  B  is 
has  a  higher  NPV  over  Case  A  as  well  as  the  potential  for  amortiza¬ 
tion  of  capital  investment  in  a  very  short  time  (1  year).  Therefore, 
Case  B  is  more  attractive  for  the  plant  operator. 

The  cumulative  cash  flow  for  each  of  these  cases  is  presented  in 
Fig.  16. 

4.2.2.  Parametric  study  and  sensitivity  analysis  to  variations  in  key 
parameters 

As  some  critical  parameters,  such  as  the  price  of  biomass  and 
the  cost  of  purchasing  C02  allowances,  exhibit  variations  based 
on  market  conditions  and/or  high  uncertainty,  a  sensitivity  analy¬ 
sis  should  be  performed  in  order  to  evaluate  the  impact  of  any 


Fig.  16.  Cumulative  cash  flows,  expressed  as  NPV  per  project  year. 
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(a)  (b) 


Fig.  17.  Sensitivity  analysis  of  NPV  to  changes  in  biomass  cost  for  (a)  Case  A  (pellets)  and  (b)  Case  B  (baled  and  shredded  material). 


C02  emission  allowance  cost  (€/t) 

(a) 


(b) 


Fig.  18.  Sensitivity  analysis  of  NPV  to  changes  in  the  cost  of  the  purchase  of  C02  emission  allowances  for  (a)  Case  A  and  (b)  Case  B. 


changes  to  their  vales  as  to  the  co-firing  project  profitability,  as  ex¬ 
pressed  by  the  NPV.  Therefore,  the  variation  in  NPV  of  investment 
is  calculated  in  relation  to  the  change  in  the  price  of  biomass  (Case 
A,  Fig.  17a  and  for  Case  B,  Fig.  17b)  and  the  change  in  the  cost  of 
purchasing  C02  allowances  (Case  A,  Fig.  18a  for  Case  B,  Fig.  18b). 

From  the  above  it  seems  that  for  Case  A,  a  fuel  cost  over  173  €/t, 
makes  the  investment  unprofitable.  On  the  other  hand,  Case  B  is 
profitable  for  all  investigated  ranges  of  fuel  prices.  Moreover,  lower 
costs  of  C02  emission  allowances  does  not  seem  to  threaten  the 
viability  of  investment,  while  the  expected  increase  in  the  coming 
years  yields  higher  NPVs  and  makes  co-firing  even  more  profitable 
for  both  cases. 

5.  Conclusions 

In  this  paper,  the  numerical  results  for  a  parametric  investiga¬ 
tion  for  the  operation  of  a  Greek  lignite-fired  power  plant  which 
co-fires  biomass  at  a  thermal  share  of  10%  are  presented  in  detail, 
while  a  first  attempt  to  compare  the  financial  viability  of  two  bio¬ 
mass  handling  and  feeding  scenarios  is  performed.  Overall,  the 
most  important  operational  parameters  of  the  boiler  are  unaffected 
by  this  co-firing  concept  (low  thermal  loading),  despite  slightly  in¬ 
creased  furnace  exit  temperatures  for  some  investigates  scenarios; 
the  observed  increases  are  not  so  significant  for  one  to  expect 
intensifying  slagging/fouling  phenomena  related  to  the  lignite  par¬ 
ticles  in  comparison  with  the  reference  case  of  operation.  The  total 
heat  transfer  to  the  furnace  wall  is  likewise  unaffected.  CFD  analy¬ 
sis  suggests  that  a  potential  benefit  of  co-firing  conditions  is  the 
decrease  of  NOx  emissions  up  to  10%,  due  mostly  to  the  lower 


nitrogen  content  of  the  biomass  fuel  and  the  respective  mechanism 
for  fuel  NO*  formation.  From  this  parametric  study,  it  is  deduced 
that  in  order  to  achieve  satisfactory  combustion  in  the  case  of  bio¬ 
mass  and  avoid  increased  unburnt  losses,  a  grinding  of  the  biomass 
fuel  in  mean  diameters  of  0.5-1. 0  mm  is  required;  for  larger  bio¬ 
mass  particles  (3  mm),  a  high  portion  of  unburnt  char  in  the  ashes 
is  expected.  Additionally,  large  biomass  particles  result  in  high 
temperatures  at  boiler  heights  lower  than  12  m,  which  may  cause 
issues  to  the  boiler  operation,  especially  if  the  substitution  per¬ 
centage  is  increased  to  values  higher  than  10%.  In  term  of  injection 
points,  the  most  efficient  case  is  to  inject  the  biomass  particles 
from  the  vapour  burners,  though  the  upper  main  burners’  level  is 
another  possible  alternative. 

Taking  into  account  possible  market  prices  for  fuels  and  for  the 
equipment  required  for  the  co-firing  retrofit,  an  economic  evalua¬ 
tion  of  the  different  cases  suggests  that  Case  B,  although  more 
technically  complicated  and  with  a  higher  investment  cost,  is  more 
profitable.  This  is  mostly  due  to  the  lower  fuel  cost  of  baled  and/or 
shredded  material  compared  to  pellets.  The  effect  of  variations  of 
several  key  parameters  such  as  the  price  of  biomass  and  the  cost 
of  purchasing  C02  allowances  in  the  viability  of  the  investment 
also  supports  this  claim. 

The  results  of  the  economic  analysis  suggest  that  an  additional 
investment  for  dedicated  biomass  burners  is  expected  to  result  in 
profitable  projects,  while  leading  to  further  optimisation  of  the  com¬ 
bustion.  Based  on  the  results  of  this  investigation,  a  suitable  level  for 
biomass  injection  lies  at  heights  higher  than  the  upper  main  burner 
level  ( 1 6.5  m),  so  as  to  ensure  reduced  unburnt  losses  in  the  bottom 
ash  and  satisfactory  burnout  of  the  biomass  at  the  outlet  of  furnace. 
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